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ABSTRACT 

I Zw 18 is the most metal poor star-forming galaxy known and is an ideal laboratory 
to probe stellar evolution theory at low metallicities. Using archival HST WFPC2 imag- 
ing and FOS spectroscopy we were able to improve previous studies. We constructed a 
continuum free Hell A4686 map, which was used to identify Wolf-Rayet (WR) stars re- 
cently found by ground-based spectroscopy and to locate diffuse nebular emission. Most 
of the Hell A4686 emission is associated with the NW stellar cluster, clearly displaced 
from the surrounding shell-like [O III] and Ha emission. We found evidence for He II 
sources, compatible with 5-9 WNL stars and/or compact nebular Hell A4686 emission, 
as well as residual diffuse emission. Only one of them is outside the NW cluster. We 
have done an extensive comparison between our results and the recent ground-based 
data used by Izotov et al. (1997) and Legrand et al. (1997) to identify WN and WC 
stars in I Zw 18. The differences between the various data may be understood in terms 
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of varying slit locations, continuum fits, and contamination by nebular lines. We have 
calculated evolutionary tracks for massive stars and synthesis models at the appropri- 
ate metalhcity {Z ~ 0.02 Zq). These single star models predict a mass limit Mwr ~ 
90 Mq for WR stars to become WN and WC/WO. For an instantaneous burst model 
with a Salpeter IMF extending up to Myp ~ 120-150 Mq our model predictions are in 
reasonable agreement with the observed equivalent widths. Our model is also able to 
fully reproduce the observed equivalent widths of nebular Hell A4686 emission due to 
the presence of WC/WO stars. This quantitative agreement and the spatial correlation 
of nebular Hell A4686 with the stellar cluster and the position of WR stars shown from 
the ground-based spectra further supports the hypothesis that WR stars are responsible 
for nebular He II emission in extra-galactic H II regions. 

Subject headings: galaxies: individual (I Zw 18) — galaxies: stellar content — stars: 
Wolf-Rayet 
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1. Introduction 

I Zw 18 (= Mrk 116 = UGCA 166) is described 
in Zwicky's catalog of compact galaxies (Zwicky 
1971) as a pair of interconnected blue compact, 
presumably young galaxies. This otherwise un- 
remarkable galaxy attracted immediate attention 
when Searle &; Sargent (1972) established its ex- 
tremely low oxygen abundance of only 2% the 
solar value. I Zw 18 became one of the proto- 
types of the "flashing galaxies" (Searle, Sargent, 
&; Bagnuolo 1973), whose blue colors suggest a 
brief period of intense star formation. Subse- 
quently this galaxy class was generally referred 
to as "star burst galaxies" (Weedman 1987), with 
I Zw 18 belonging to the sub-category of blue 
compact dwarfs (Thuan 1991). A review of its 
most pertinent parameters was given by Du- 
four et al. (1996): It has a radial velocity of 
740 km s~^, is at a distance of about 10 Mpc, 
and has an extinction of E{B — V) = 0.07. 

I Zw 18's extreme metal deficiency was con- 
firmed in several follow-up studies (e.g., Kunth &; 
Sargent 1986; Campbell 1990; Pagel et al. 1992; 
Kunth et al. 1994; Stasihska & Leitherer 1996). 
Although strong efforts to detect other galax- 
ies with very low metal content have been un- 
dertaken (e.g. Terlevich, Skillman, & Terlevich 
1995), I Zw 18 is still the record holder. Its low 
metallicity is of particular importance for stud- 
ies of the primordial helium abundance (Izotov, 
Thuan, Sz Lipovetsky 1996; Olive, Skillman, & 
Steigman 1997; Skillman, Terlevich, & Terlevich 
1998), and for properties of chemically unevolved 
galaxies in general. Despite the low oxygen abun- 
dance, however, there is clear evidence for stellar 
nucleosynthetic processing from its carbon abun- 
dance so that I Zw 18 is not a primordial galaxy 
(Garnett et al. 1997). The same suggestion was 
made earlier by Thuan (1983) on the basis of 
surface photometry of its underlying stellar pop- 
ulation. 

I Zw 18 is an ideal laboratory to probe stel- 
lar evolution theory in a low-metallicity environ- 



ment. The spatial resolution of WFPC2 on the 
Hubble Space Telescope (HST) allows the detec- 
tion of individual stars down to the ~10 Mq 
range (Hunter &; Thronson 1995; Dufour et al. 
1996). Analysis of deep WFPC2 images clearly 
demonstrates ongoing massive-star formation. Ad- 
ditional, indirect evidence for a population of 
massive stars comes from bubbles and outflows 
of the interstellar medium, indicative of powerful 
hot-star winds and supernovae (Martin 1996). 

Wolf-Rayet (WR) stars are the evolved, less 
massive descendants of previously massive O stars 
(Maeder &; Conti 1994). Stellar evolution the- 
ory predicts few WR stars to form in a low- 
metallicity environment such as in I Zw 18 (e.g., 
Meynet 1995). Therefore recent reports (Izotov 
et al. 1997, hereafter IFGGT; and Legrand et 
al. 1997, hereafter LKRMW) of WR star detec- 
tions in I Zw 18 have far-reaching implications for 
stellar evolution theory. Furthermore I Zw 18 is 
probably the best studied object among approxi- 
mately 50 to 60 extra-galactic H II regions show- 
ing nebular Hell A4686 emission (Campbell et al. 
1986; Izotov, Thuan & Lipovetsky 1994, 1997) 
whose origin is still poorly known and subject 
to recent debates. Various mechanisms reach- 
ing from photoionization by hot massive stars, 
photoionization by X-ray binaries to shocks have 
been put forward to explain the high excitation 
(Garnett et al. 1991). Understanding the origin 
of nebular Hell A4686 in low-metallicity galax- 
ies has implications for our knowledge of the 
far-UV flux of galaxies and may prompt a re- 
examination of the importance of young galaxies 
to the ionization of QSO absorption line systems 
and the ionization of the intergalactic medium 
(cf. Garnett et al. 1991; Pettini et al. 1997). 

WR stars can be recognized, e.g., via their 
broad (~100 A) emission bump around A4686, 
which is generally a blend of H, He II and sev- 
eral metal lines. I Zw 18 has been known to 
show narrow, nebular Hell A4686 emission but 
the broad WR emission bump was not known be- 
fore the findings of IFGGT and LKRMW. The 
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blending of several stellar and nebular emission 
lines around 4686 A and the complex spatial mor- 
phology in ground-based data make it challeng- 
ing to disentangle stars and gas and to derive the 
WR content and the nebular properties. This is 
evident in the apparently discrepant results of 
IFGGT and LKRMW. The superior spatial res- 
olution of HST may prove useful in resolving the 
discrepancy. 

In this paper we report on an HST archival 
project to study the spatial morphology of the 
highly ionized gas in I Zw 18, as observed pri- 
marily in Hell A4686. Our goals are to recon- 
cile existing ground-based spectroscopic observa- 
tions, to test stellar evolution models for WR 
stars in an extreme chemical environment, and 
to elucidate the nature of nebular He II emis- 
sion in extragalactic H II regions. The paper is 
organized as follows. In §2 and §3 we describe 
the data used in this work and their analysis. In 
§4 we discuss the recent spectroscopic detection 
of Wolf-Rayet stars in connection with the HST 
imagery. The interpretation and a comparison 
with Wolf-Rayet star models are presented in §5. 
Finally, §6 summarizes our conclusions. 

2. Data and reduction 
2.1. Observation log 

We have retrieved HST archival data of I Zw 18 
obtained during Cycles 4 and 6. The data con- 
sist of both Wide Field/Planetary Camera 2 
(WFPC2) images and Faint Object Spectrograph 
(FOS) spectra. The existence of such an exten- 
sive collection of data including images and spec- 
tra in several passbands make it possible to im- 
prove previous analysis, which were done using 
the available data at that time and which inter- 
preted each data set individually (e.g. Hunter 
& Thronson 1995; Dufour et al. 1996). In Ta- 
ble 1^ we list the image root names (column 1), 
proposal identification number (column 2), the 
epoch of observation (column 3), the configura- 
tion (column 4), and the optical element (col- 



umn 5) for each set of data retrieved from the 
HST archive. In column 6 we give the spectral 
range of each optical element. The entry is the 
full width at half maximum (FWHM) for filters 
and the free spectral range for gratings. In the 
following subsections we describe several meth- 
ods of data reduction which were tested using 
several different images and spectra. 

2.2. Pipeline processing 

All data were processed through the stan- 
dard pipeline calibration. Calibration of the FOS 
spectra in the Routine Science Data Process- 
ing (RSDP) system most notably includes detec- 
tor background and scattered light subtraction, 
Hatfield corrections, computation of wavelengths, 
and conversion from count rates to fluxes. The 
WFPC2 data were processed through the pipeline 
where bias, dark, and flatfielding corrections were 
performed and photometry keywords were calcu- 
lated. 

Since the expected flux levels are close to the 
background and the noise level of the detector, 
additional data reduction and analysis proce- 
dures on top of the standard pipelines processing 
were required. These are described in the follow- 
ing subsections. All steps were performed using 
standard tasks in IRAF/STSDAS. 

2.3. Cosmic-ray and hot pixel removal 

Multiple images were simultaneously co-added 
and cosmic-rays removed using the IRAF task cr- 
rej. We set the algorithm to estimate the initial 
value to 'median' in order to eliminate the pos- 
sibility of inadvertent rejection of true counts. 
Our study approaches the limits of HST detec- 
tion, and high cosmic-ray flux values would bias 
the mean towards a higher number. Likewise by 
setting the a threshold higher, we ensured the ac- 
curacy of the cosmic-ray rejection task. Surface 
plots of all supposed Hell A4686 emission were 
visually inspected for the characteristic profile 
of a cosmic-ray hit. The probability that three 



4 



cosmic-rays can hit the same pixel in three dif- 
ferent images was tested by using dark images 
taken at the same epoch and with similar ex- 
posure time. The dark images were cleaned us- 
ing the same procedure as described above. We 
found that only 6 pixels remained contaminated 
by cosmic-rays (total number of pixels in each im- 
age = 6.4 X 10^). Considering that this analysis 
will be critical when we discuss the Hell A4686 
sources in §3, it is important to notice that this 
number is even lower if we take into account that 
the area of the CCD where the helium sources are 
detected is much smaller (total number of pixels 
= 4.7 X 10^). Therefore, our technique of elimi- 
nating cosmic-rays proved to be very efficient and 
adequate for this study. 

All pixels with apparent Hell A4686 emission 
were also individually inspected and compared 
with the known coordinates of hot pixels on the 
chip to eliminate any erroneous detections. 

2.4. Background subtraction 

Background subtraction was also performed 
on each image. Automatic background subtrac- 
tion through crrej which simply subtracts the 
mode of the image was not used because we de- 
termined that this method generally subtracted 
values much larger than the true background. In- 
stead, the background value was determined by 
computing the average of the counts calculated 
from the median of 13 boxes of 100 square pixels 
placed in relatively blank regions of each image. 
This procedure was performed on all images and 
the corresponding background value subtracted. 

2.5. Flux calibration 

All images were then converted into absolute 
flux units. Each image was multiplied by the ap- 
propriate PHOTFLAM value (erg s"^ cm^^ A"^) 
which is defined to be the mean flux density that 
produces a count rate of 1 per second (DN) with 
the HST observing mode used in each observa- 
tion. The final flux calibrated image was created 
by dividing by the exposure time. 



The reliability of this process was exhaus- 
tively confirmed through comparisons with sev- 
eral sources of spectral data. Photometry was 
performed on the WFPC2 images and compared 
with the flux values obtained from the FOS spec- 
tra. Measurements of spectra from proposal 
6536 (y39a0304t) indicated hue fluxes of 9.50 x 
10~^^ erg s~^ cm~^ and 7.59 x 10~^^ erg s~^ cm^^ 
for Ha and [OIII] A5007 ([OIII] hereafter), re- 
spectively. These values refer to a circular en- 
trance aperture with a diameter of 0.86". We 
placed a synthetic aperture on the Ha and [OHI] 
WFPC2 images and measured the encircled fluxes. 
Photometry was then performed on the Ha and 
the continuum subtracted [OIII] WFPC2 im- 
ages generated from observations of proposal 
5434. (See below for the adopted continuum 
subtraction technique.) We measured 9.83 x 
10~^^ erg s~^ cm~^ for Ha and 6.54x 10~^^ erg cm~^ 
for [OIII] . We consider these values in good agree- 
ment, given the inherent uncertainties, most no- 
tably the unknown absolute pointing of HST. 
Typical pointing uncertainties are ±0.5". Pho- 
tometry was also performed at offset positions 
approximately 1" away from the assumed spec- 
tra locations to test the effect of HST pointing 
uncertainties. Flux differences on the order of 
50% were found. 

Total flux levels of Ha and [OIII] for the 
entire galaxy as measured from the flux cali- 
brated WFPC2 images were likewise compared 
with published values from ground-based spec- 
tra. In Table ^ we compare our flux measure- 
ments with those of several different authors. 
Our WFPC2 fluxes are for the entire galaxy. 
The published ground-based measurements have 
a spread of nearly an order of magnitude, with 
the extreme values bracketing the WFPC2 val- 
ues. Most likely, this reflects uncertainties in 
the calibration of ground-based data, different 
filter systems, and the difficulty to measure faint 
extended emission on ground-based data. Our 
values are closer to Dufour & Hester (1990) val- 
ues which are based on imaging (see also Martin 
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1996). 

Wc arc mostly concerned about the inter- 
nal uncertainty of the WFPC2 photometry since 
we combine different datasets to construct pure, 
continuum-free line images. To address this con- 
cern, flux levels of multiple F555W images were 
compared. We used a WF3 and a PC image from 
proposal 5434 and a PC image from proposal 
5309. Fluxes at discrete locations as well as the 
total flux of the main body of I Zw 18 were com- 
pared in this filter. The total fluxes from these 
three images are 1.11 x 10^^^ erg s^^ cm~^ 
1.14 X 10~^^ erg s~^ cm~^ A~^, and 1.10 x 
10~^^ erg s~^ cm~^ A~^. We conclude from these 
cross-checks that the internal WFPC2 flux cali- 
bration is accurate to about 3%. 

2.6. Continuum subtraction 

Pure emission-line images were produced by 
subtracting the underlying continuum from the 
appropriate narrow-band image. A pure oxygen 
map was created by subtracting a scaled con- 
tinuum image compatible with the wavelength 
range of the F502N filter. We also created an 
Ha map. The Ha emission was taken as the flux 
in the F658N filter as the Ha line is redshifted to 
'--'6580 A, which is within the bandwidth of the 
F658N filter (FWHM = 28.5 A). A very conser- 
vative upper limit of 10% for the continuum con- 
tribution was estimated from FOS spectra. [Nil] 
A6548/83 is outside the passband of the F658N 
filter. Hence we use the flux calibrated F658N 
image as an Ha map without any further correc- 
tion. 

A continuum-free spatial map of Hell A4686 

was created using the F336W image since this 
image contains no significant line emission. Based 
on the FOS spectra of proposal 6536, the ratio of 
the continuum fluxes of the F469N and F336W 
fllters at their effective wavelengths is 0.35. The 
F336W image was scaled by the appropriate fac- 
tor and then subtracted from the F469N for a 
pure Hell A4686 map. It is important to notice 
that the FOS spectra were centered on a small 



region (diameter = 0.86") on the southwest side 
of the NW region of I Zw 18 and not necessarily 
represent the entire NW region. In order to check 
this we have compared FOS fluxes with ground 
based fluxes (F. Legrand 1997, private commu- 
nication). We found that the FOS spectra are 
quite representative of the entire NW region, i.e. 
larger scales. However, in smaller scales the same 
cannot be assured. To put the following discTis- 
sion into perspective, we note that the pure Hell 
A4686 flux is about 60% of the continuum flux 
within the F469N filter. 

The emission-line images used for the anal- 
ysis were obtained as described in the previous 
paragraphs. The adopted continuum subtraction 
technique was the result of extensive tests using 
diflFerent methods. Ideally, one would like to use 
a narrow-band image in the adjacent, line- free 
continuum for the continuum subtraction. Given 
HST's filter set, such data do not exist. Errors 
in the continuum subtraction process are intro- 
duced by using (i) continuum images which are 
contaminated by lines and (ii) continuum images 
whose effective wavelength is very different from 
that of the line image. In the following we dis- 
cuss our experiments to evaluate the trade-offs 
between (i) and (ii). Since our main interest is in 
the He II map, we restrict our discussion to the 
continuum correction of the F469N filter. 

Continuum maps were created by correcting 
the F555W images for emission lines. As deter- 
mined through equivalent width measurements of 
FOS spectra from proposal 6536 (y2f90402t), the 
emission lines in this filter contribute about 40% 
to the total flux of the F555W fllter, and so this 
fraction was subtracted to generate the contin- 
uum map. These corresponding continuum im- 
ages were then subtracted from the F469N image 
to create a He II map. However, the large correc- 
tion required to account for the line contamina- 
tion in the F555W filter and its probable spatial 
dependence introduce a large uncertainty. The 
resulting He II maps turned out to be unreliable. 

Other experiments involved various combina- 
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tions of the F450W and F439W images together 
with F336W and F555W to inter/extrapolate the 
continuum level at 4686 A. Again, we noticed 
the uncertain line contamination correction for 
the continuum filters. This was tested by de- 
termining the location of the continuum using 
F336W (no significant contamination), F439W 
and F450W (some contamination), and F555W 
(strong contamination). Our tests showed that 
by far the more important error source in cre- 
ating a He II map is the line contamination in 
the continuum filters and their spatial variation. 
In contrast, the stellar continuum slope is not 
strongly variable. This is immediately obvious: 
stars contributing to the optical/near-ultraviolet 
light have spectral types O to early A. Their 
(B — V) colors are in the narrow range of -0.3 to 
-|-0.1 so that the continuum slope is fairly con- 
stant and has little spatial dependence. 

We conclude from these simulations that our 
Hell A4686 map is well-suited to detect any ex- 
cess He II emission but that its absolute flux cal- 
ibration must be considered with care. We adopt 
a conservative error of 30%. Next we discuss the 
Ha, [OIII], and Hell A4686 maps. 

3. Data Analysis 

In Fig. I we show the WFPC2 V (F555W), 
Ha (F658N) and (F502N) images of I Zw 18The 
total Ha and [OIII] fluxes are given in Table |, 
which was described in §2.5. The Ha and [OIII] 
maps have a very similar morphology with coin- 
cident peaks of emission and similar filamentary 
structure. The ionized gas distribution is dis- 
played in forms of 'bubbles' probably driven by 
stellar winds and super novae. The nebular peaks 
of Ha and [OIII] are spatially offset from the U, 
V and I continuum. We refer to Hunter & Thron- 
son (1995) and Martin (1996) for a more detailed 
analysis of the ionized gas and its relationship to 
the stars. 

The spatial map of continuum- free Hell A4686 
was carefully analyzed and 25 individual pixels of 



He II emission greater than 3cr above the 'back- 
ground' level have been identified within an area 
of 9.9 arcsec^ (62 pixels > 2a). This area is cen- 
tered on the NW component of the galaxy and 
is presented in Fig. |2| together with the V im- 
age. The darkest pixels in Fig. |2| have fiuxes 
higher than 3a. We have verified that these He II 
sources are neither contaminated by cosmic-rays 
nor by hot pixels (see §2.3 for more details on 
cosmic-ray removal). As a final check we identi- 
fied each Hell A4686 emission in the three indi- 
vidual images used in crrej. Since the probability 
that three cosmic-rays hit the same pixel in three 
different images is very low (see §2.3) this proce- 
dure guarantees that each emission detected in 
the three images are not cosmic-rays, i.e. pixels 
that were not detected in each of the three images 
were considered contaminated by cosmic-ray. 

From the present data we have no direct means 
of distinguishing nebular from stellar Hell A4686 
emission. Stars of the types Of, Ofpe/WN, 
WNL, and WNE show Hell A4686 emission formed 
in the expanding stellar wind. A fraction of the 
broad blended emission of C III/C IV AA 4640- 
4650 and Hell A4686 in WC and WO stars may 
also be detected in the F469N filter. The average 
stellar He II line luminosity of WR stars of the 
WN sequence is L4686(WNE) = (5.2 ±2.7) x 10^5 
erg s-^ and L4686(WNL) = (1.6 ± 1.5) x 10^*^ 
erg (see Schaerer &: Vacca 1998 for line lu- 
minosities of other types). Placing a WNL star 
at the distance of I Zw 18 (10.8 Mpc) yields a 
typical line flux of 1.15 x 10~^^ erg s~^ cm~^. 
The He II fiher (F469N) has a width of 24.9 A 
and therefore the average line fiux of a WNL star 
detected with F469N would be expected to be 
of the order of 4.6 x 10"^® erg s~^ cm~^ A~^ 
(average line fiux over the FWHM of the fil- 
ter). However, we should take into account the 
fact that Galactic WN stars have Hell A4686 
13 A<FWHM<63 A and WN stars in the SMC 
have 17 A<FWHM<33 A (Conti k Massey 1989; 
Conti et al. 1989). Considering that I Zw 18's 
stars are more similar to SMC's stars than to 
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the ones in the Milky Way, we conclude that 
the average line flux detected with the F469N 
(FWHM=24.9 A) contains most of the true emis- 
sion. 

The total He II flux coming from helium 
sources (>3cr) in the region presented in Fig. ^ 
is 3.97 X 10~^^ erg s~^ cm~^ A.~^, equivalent to 
~ 9 WNL stars 0. However, this flux is dis- 
tributed over several regions. For instance, 12 
individual pixels are above the 5a level and have 
fluxes higher than 1.0 x 10~^® erg s~^ cm~^ A~^. 
We have identified 6 of these sources as pos- 
sible WR stars (>10(t). The pixels marked 
with WR and "WR?" in Fig. | have a flux 
of 1.3 X lO-^'^ erg s~^ cm'^ A'^ (~ 3 WNL 
stars) and 3.4 x 10"^^ erg s~^ cm~2 A^^ (~0.7 
WNL star), respectively. The pixels marked with 
"WR??" have fluxes of 1.7, 3.0, 5.3 and 9.2 
xlO~^^ erg s~^ cm~^ A~^, which is equivalent 
to 0.4, 0.7, 1.2 and 2 WNL stars (4.3 WNL in 
total), respectively. However, only one of these 
pixels (0.4 WNL) is classified as an unambigu- 
ous detection. The other three pixels are present 
in two of the images but are seen as faint detec- 
tions in one of the images. Therefore, we can- 
not exclude the possibility of contamination by 
cosmic-rays (see §2.3 for more details on cosmic- 
ray removal). 

Although the fluxes from the individual He II 
sources found are compatible with fluxes of WNL 
stars, we cannot exclude sources of compact neb- 
ular Hell A4686 emission or a combination of 
stellar and nebular emission. Placed at the dis- 
tance of I Zw 18 and observed with the F469N fil- 
ter the 3 nebular Hell A4686 sources DRl, SMC 
N76, and LMC N44C studied by Garnett et al. 
(1991) have fluxes between 1.6 x IQ-^^ (N76) and 
1.8 X 10"^^ erg s"^ cm"^ A"^ (N44C), compara- 
ble to 0.4-3.6 WNL stars. Their angular diame- 

■^The flux value is not corrected for extinction. Using the 
values of Hunter & Thronson (1995) for Aa (0.21, 0.13, 
and 0.08 for F336W, F555W, and F814W, respectively) 
and assuming that the Milky Way extinction curve is ad- 
equate for dereddening, we find A4686 = 0.18. 



ters range from ~ 0.4" (DRl) to 0.03 " (N44C) at 
the distance of I Zw 18. Therefore, our He II de- 
tections are also compatible with nebular sources 
similar to Garnett's objects provided that they 
have compactness similar to N44C. 

The integrated line flux within the 9.9 arcsec^ 
region including the helium sources is found to 
be 1.67 X 10"^'^ erg s"^ cm'^ A'^ i.e. 4.16 x 
10~^^ erg s^^ cm~^. For comparison, ground- 
based observations of the narrow Hell A4686 
emission in the NW component of I Zw 18 yield 
somewhat lower values between 5.0 x 10~^^ 
and 1.3 x 10^^^ erg s~^ cm~^ (Pagel et al. 
1992; Skillmann & Kennicutt 1993, IFGGT; Izo- 
tov &: Thuan 1998), which is likely explained 
by the smaller area covered by these observa- 
tions. In our images the difference between the 
total flux and the flux from the helium sources 
is 1.1x10"^^ erg s~^ cm~^ A~^. This resid- 
ual emission is spread over the 9.9 arcsec^ re- 
gion and is probably nebular. We have also 
checked for the possibility of contamination by 
the red stellar objects found by Hunter & Thron- 
son (1995) since these objects may contribute 
flux at 4690 A but not in F336W. In Fig. | 
we show the helium sources and all the red 
stars brighter than F555W=25.0 and redder than 
F555W-F814W=0.5 (see Hunter & Thronson 
1995 for details on magnitude calculation). We 
find that 5 of these red stars coincide with the 
position where helium sources were identified. 
However, these helium sources are very faint with 
fluxes of ~ 5 — 9 X 10~^^ erg s~^ cm~^. The total 
helium emission from these sources contributes 
less than 6% to the total flux in the entire re- 
gion. 

In order to compare the features of differ- 
ent maps with the Hell A4686 map. Wide Field 
(WFC) images of Ha (u2f90205t and u2f90206t) 
and [OIII] (u2f90203 and u2f90204) were re- 
binned into Planetary Camera (PC) resolution. 
We used the task magnify in IRAF to perform the 
image interpolation. We have used 5 V images 
(u2cg0201t, u2cg0202t, u2cg0203t, u2f90104t and 
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u2f90105t) taken with the PC and WF3 in or- 
der to check the vahdity of this technique. We 
conclude that the transformed images preserve 
enough information on the general morphology of 
the maps for a qualitative inspection. However, 
the transformed image loses information during 
the image interpolation and should be used with 
caution when comparing detailed features. 

The helium peaks are located in the NW 
cluster, whereas the [O III] and Ha emission 
has its maximum in a shell-like structure sur- 
rounding the cluster. Measured in small circu- 
lar areas around the possible WR stars one finds 
/(4686)//(Ha) ^ 0.3-0.1; a fairly constant value 
of /(5007)//(Hq) ~ 1 is obtained over the en- 
tire region of the cluster. However, these values 
should be taken with caution since they were ob- 
tained using the WF images transformed to PC 
resolution. 

Several stars are identified in the area of the 
He II peaks, including the second brightest V star 
in the list of Hunter & Thronson (1995). This 
star is only 4 pixels (1 pixel=0.046") to the west 
of the WR stars and according to those authors 
it has magnitudes F555W=22.43, F336W=21.67 
(see Hunter & Thronson 1995 for details on their 
photometry) . 

We have also investigated the presence of he- 
lium sources outside the NW cluster. We found 
that one source has flux ~1.5 WNL. It is located 
in between the NW and the SE condensations. In 
Table | we summarize the helium sources identi- 
fied in I Zw 18. 

Prom the spatial maps presented in this sec- 
tion we conclude that Hell A4686 emission is 
clearly associated with the stellar cluster and 
spatially offset from the maximum nebular Ha 
and [OIII] emission. The present data alone do 
not allow us to distinguish between nebular and 
stellar Hell A4686 emission. 



4. Comparison with ground-based spec- 
troscopy 

Two very recent papers have reported dif- 
ferent numbers of WR stars in I Zw 18 using 
ground-based spectroscopy. LKRMW estimated 
from spectra taken with the 3.6m CFH Telescope 
that 1-2 WC4 or WC5 and no WN stars are 
present in I Zw 18. In contrast, IFGGT using 
the Multi Mirror Telescope found 17±4 WNL 
and 5ib2 WC4 stars. Both spectra unambigu- 
ously show broad blue and red emission bumps 
centered at ~4645-4686 A and at 5808 A, respec- 
tively, but differ in the measured line fluxes and 
the structure of the blue emission bump, as sum- 
marized below. There may be several reasons 
why their results differ. The main reason may be 
the choice of a different slit position and orienta- 
tion. The slits were centered approximately on 
the same region (central knot of NW H II region) . 
However, the position angle used by LKRMW is 
+A5° compared to -41° by IFGGT. Therefore, 
their spectra cover partially different regions of 
the galaxy. 

For the following discussion and for compar- 
isons with theoretical predictions of the massive 
star population in I Zw 18 (see §5) we com- 
plement both sets of measurements of LKRMW 
and IFGGT with equivalent width measurements 
kindly provided by P. Legrand (1997, private 
communication), and with additional H/3 equiv- 
alent widths from the literature. The data are 
summarized in Table ^. Before proceeding we 
give the following cautionary remarks: (i) LKRMW 
detect a broad component centered at ~4645 A 
(cf. "broad 4645" in Table §) and no broad emis- 
sion at A ~ 4686 A. A small change in the fit of 
the continuum level may, however, accommodate 
the existence of a broad component, (ii) IFGGT 
detect a broad "WR bump" extending over A ~ 
4619 - 4740 A ("broad blue bump" in Table |). 
As mentioned by these authors this feature can 
generally be composed of several broad WR emis- 
sion lines as well as nebular lines. The presence 
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of [Felll] A4658, [ArIV]+HeI AA4711,4713, and 
[ArlV] A4740 can be suspected in their spectrum 
(the latter two hnes are present in LKRMW). 
When present, these nebular emission lines typi- 
cally have line intensities of 0.01 to 0.03 H/3 in 
total, as can be seen from the data of Izotov, 
Thuan k Lipovetsky (1994, 1997). The contri- 
bution due to WR stars may therefore be overes- 
timated by up to a factor of two. (iii) Both the 
intensity of CIV A5808/H/3 and the CIV A5808 
line flux integrated over the respective regions 
differ approximately by a factor of two between 
LKRMW and IFGGT. This is very likely ex- 
plained by the complex spatial distribution of the 
II/5 emission as seen from the HST maps and the 
different adopted slit positions/orientations, (iv) 
The observations of the nebular Hell A4686/H/3 
intensity (~0.040) of LKRMW and IFGGT are 
in agreement but are slightly larger than previ- 
ous measurements (Hell A4686/H/3 0.032; e.g. 
Pagel et al. 1992; Skillman & Kennicutt 1993). 
(v) The H/3 equivalent widths measured for the 
NW HII region differ considerably between var- 
ious authors. Values between 56 and 127 A are 
found by Pagel et al. (1992), Skillman & Kenni- 
cutt (1993), LKRMW and IFGGT. 

Comparing the ground-based results of LKRMW 
and IFGGT with our results (see Table ^ for a 
summary our final numbers) it is possible that 
LKRMW have underestimated the number of 
WN stars while IFGGT et al. have overestimated 
it. The HST helium filter used, F469N, covers 
the region of 4683-4707 A and therefore the he- 
lium emission line which is redshifted to 4698 A 
is well centered in the filter passband. There is 
no contamination by the lines that are blended 
in the 4645 A bump; the two closest lines, [Felll] 
A4658, [ArlV] A4711 are redshifted to 4669 A and 
4722 A, respectively, lying outside the helium fil- 
ter. Therefore the helium emission detected can 
be either stellar or nebular in origin. However, 
we cannot exclude the presence of WC stars since 
there are no HST data available for the region of 
the 5808 A bump. 



The absence of WN stars was inferred by 
LKRMW from the absence of broad 4686 A emis- 
sion underlying the narrow 4686 A. However, 
a small change in their fit of continuum level 
would provide enough flux to account for a few 
WN stars. It is also possible that no WN stars 
were detected because none were in the region 
where their spectra were taken; their slit width 
is 1.52" and could be centered just below the re- 
gion where the stars are (as discussed at the end 
of this section). 

The number of WNL stars derived by IFGGT 
is most likely overestimated for the following rea- 
sons: (i) The contribution of WC stars to the 
blue WR bump (A ~ 4686 A) has not been ac- 
counted for, and (ii) the broad Hell A4686 bump 
may be contaminated by other nebular emis- 
sion lines as discussed above. Indeed the ob- 
served flux ratio /(4686)//(5808) ^ 2 is typical 
for WC4-5 stars (cf. Smith et al. 1990; Schaerer 
& Vacca 1998) and the contribution of CIII/IV 
A4650 from the WC stars to the 4686 A bump 
should be subtracted to determine the remain- 
ing luminosity from WNL stars. Using the re- 
cent compilation of WR line luminosities from 
Schaerer & Vacca (1998) we obtain: A'wc4 ~ 4, 
and AVnl = 4. Note that the main difference 
between the value A'wnl = 17 ± 4 obtained by 
IFGGT is due to (i). If (ii) were also taken into 
account, the inferred number of WNL stars may 
be even lower. In any case it appears that the 
differences between the estimates of the number 
of WR stars derived by IFGGT and LKRMW 
are not significant if uncertainties are taken into 
account. Future observations with higher dispe- 
rion and signal-to-noise will be necessary to ad- 
dress these issues with more accuracy. In abso- 
lute terms, the number of WNL stars in partic- 
ular may be more uncertain given the relatively 
large spread of Hell A4686 line luminosities of 
individual WNL stars (Schaerer & Vacca 1998, 
their Fig. 1). 

LKRMW and Izotov & Thuan (1998) have an- 
alyzed the spatial location of several emission fea- 
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tures in I Zw 18. LKRMW have found that the 
bumps at 4645 A and 5820 A peak between 1" 
and 2" SW from the central cluster. They have 
also found a spatial correlation between WC fea- 
tures and nebular Hell A4686. In order to re- 
produce their results we have added (IRAF task 
improject) the flux of the HST images within 1.5" 
(their slit width). In Fig. |^ we present this re- 
sult for the helium and F555W images (pixel 
size is 0.046"). Three helium peaks are eas- 
ily identified; two of them are the "WR" and 
"WR?" and the third peak is the one marked 
with "WR??" in Fig. |2|. Comparing the helium 
peaks with the F555W peaks we conclude that 
the "WR??" peak is probably the one claimed 
by LKRMW as the peak located SW of the cen- 
tral cluster. Due to the lack of resolution present 
in LKRMW spectra, the central cluster is prob- 
ably centered in between the three peaks in the 
F555W (around pixel 50). "WR??" is 1" away 
from pixel 50, which agrees with the distance 
found by LKRMW. However, as we discussed 
in §3 there is the possibility of cosmic-ray con- 
tamination at the "WR??" position. Our re- 
sults suggest that LKRMW's slit location prob- 
ably missed the "WR" and "WR?" peaks, oth- 
erwise one bump should have been detected on 
the NE of the central cluster. The two peaks, 
"WR" and "WR?" , cannot be resolved with the 
ground-based telescopes since they are only 0.1" 
apart situated on the NE of the central cluster. 
The low spatial resolution and uncertainties in 
pointing of the ground-based spectra limit the 
comparison with HST data. 

The shift between the hydrogen emission and 
the continuum found by LKRMW agrees well 
with the HST images. Our spatial maps are also 
compatible with the spatial distribution found by 
Izotov k Thuan (1998). The results of IFGGT 
do not provide information about the position of 
the WR features. 



5. Model comparison and interpretation 

In order to compare all the data with appro- 
priate model predictions of massive star popula- 
tion at the low metallicity of I Zw 18 we have first 
calculated a set of new stellar evolution tracks 
and then performed evolutionary synthesis cal- 
culations using these tracks. Preliminary results 
have been discussed in Schaerer (1997). The in- 
gredients of the stellar evolution models and a 
discussion of the results from these calculations 
are presented in the next subsections. 

5.1. Stellar evolution models 

Evolutionary tracks for stars between 25 and 
150 Mq (Mini = 25,40,60,85,120,150) at the 
metallicity Z = 0.0004 ~ 1/50 Zq have been cal- 
culated with the Geneva stellar evolution code 
adopting the same ingredients as Meynet et al. 
(1994). These calculations extend the previously 
available Geneva sets (0.001 < ^ < 0.1) down to 
the metallicity of I Zw 18. The adopted mass-loss 
rate and its metallicity dependence are the main 
factors determining the evolution of the most 
massive stars. As shown by Maeder &: Meynet 
(1994) the high mass-loss rates adopted in the 
models of Meynet et al. (1994) reproduces many 
observational properties of individual WR stars 
and O star populations at different metallicities, 
although direct measurements of WR mass-loss 
rates at Zq suggest lower values (Leitherer et 
al. 1997). We have therefore adopted the same 
mass-loss prescription and scaling with metallic- 
ity. In addition we have also computed evolu- 
tionary tracks using mass-loss rates for OB stars 
determined from the recent results of Lamers & 
Cassinelli (1996), which relies on the observed 
wind-momentum-radius relation by Kudritzki et 
al. (1995). This prescription also explicitly in- 
cludes a metallicity dependence, based on obser- 
vations of O stars in the Galaxy, the LMC, and 
the SMC. 

Qualitatively, the evolutionary tracks at Z = 
0.0004 reproduce the properties found earlier (cf. 
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Meynet et al. 1994). Due to the diminished mass 
loss only the most massive stars are predicted to 
evolve to WR stars. The (initial) mass limit of 
single WR stars is found to be Mwr ~ 90 Mq, 
compared to 61 Mq at Z = 0.001 (Maeder k 
Meynet 1994). Both the 120 and 150 M© mod- 
els evolve through the WN phase prior to be- 
coming WC/WO stars shortly before the end of 
He-burning. At the entry to this phase the ra- 
tio of the surface abundances is (C-|-0)/He > 1 
by number since the core is revealed in the late 
stage of He-burning. Following Smith & Maeder 
(1991) these stars would be classified as WO. For 
stars with initial masses ^ 80 Mq the use of 
the mass- loss prescription of Lamers & Cassinelli 
(1996) leads to essentially identical mass loss dur- 
ing the main-sequence (MS) evolution as with 
the adopted high mass-loss rates. This reflects 
the weak metallicity dependence of mass loss at 
high luminosity (see Lamers & Cassinelli 1996). 
The subsequent evolution, the WR mass limit, 
and other relevant properties are therefore nearly 
identical to the high mass loss case, although 
they are of course still subject to uncertainties 
in post-MS mass-loss rates. In summary, our 
models predict WNL stars, a short WNE phase, 
and highly evolved WC /WO stars for single stars 
with initial masses Mini ^ 90 Mq at the metal- 
licity of I Zw 18. 

5.2. Evolutionary synthesis models — com- 
parisons with observations 

Subsequently we have calculated evolution- 
ary synthesis models using the stellar tracks de- 
scribed above. The calculations are done with 
the models of Schaerer & Vacca (1998) where a 
detailed description of the input physics can be 
found. The most relevant features are the use 
of spherically expanding non-LTE atmospheres 
to describe the ionizing fluxes and the synthesis 
of WR emission features in the optical spectrum 
based on the recent compilation of line luminosi- 
ties by Schaerer & Vacca (1998). We consider 
a power-law IMF (e.g. Salpeter) with a variable 



upper mass cut-off Mup, and instantaneous star- 
formation at time t = ("instantaneous burst"). 
If the evolution of the observed WR stars in 
I Zw 18 is described by the new single star mod- 
els, Mup > Mwr ~ 90 Mq is required. Such a 
high value is very well compatible with observa- 
tions of massive stars in the Magellanic Clouds 
(Puis et al. 1996). The choice of the lower mass 
cut-off does not affect our results. 

In Fig. I we show the predicted WR/ (WR-l-0), 
WNL/(WR+0) and WC+WO/(WR+0) num- 
ber ratios for M^p =150 Mq as a function of the 
age of the burst. WR stars appear after ~2.5 Myr 
and their presence lasts approximately 1 Myr. 
The maximum WR/(WR-F0) ratio is ~0.02. 
As expected from models at Z = 0.001, WN 
stars dominate the WR population; the predicted 
WC/WN ratio is ^ 0.15. The predicted WR hues 
(line intensities with respect to H/3 and equiv- 
alent widths) and the prediction for the nebular 
Hell A4686 emission are plotted in Fig. ^ as a 
function of time. Also shown are the observed 
quantities from IFGGT and LKRMW summa- 
rized in Table It is evident that the spatial 
complexity of the gas and the displacement be- 
tween stars and gas cause strong variations in the 
stellar and nebular emission strengths. This also 
explains the large variations of W(H(3) shown in 
Table ^ On the other hand, the WR features 
and, interestingly, also nebular Hell A4686, fol- 
low fairly closely the spatial distribution of the 
continuum emission (see IFGGT, LKMRW, and 
the HST images). Therefore, we prefer to use 
the equivalent widths of the WR features as the 
physically most meaningful quantity for quanti- 
tative comparisons. Taking the largest W(H.f3) 
value (127 A) from the literature, a lower limit 
of ~3 Myr can be estimated for the age of the 
NW region from our instantaneous burst mod- 
els. 

Fig. ^b shows that the observed line intensi- 
ties exceed all predicted values by a large factor. 
This is expected due to the complex morphology 
of the interstellar gas. In contrast, the predicted 
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equivalent widths fall in the range of the observed 
values from LKRMW. We have plotted the pre- 
diction for the CIII/IV A4650 + He II A4686 
bump using the calibration from Smith (1991) 
(curves with larger values in Fig. and the 
sum of NIII A4640 and CIII/IV A4650 from the 
compilation of Schaerer & Vacca (1998) (curves 
with lower values in Fig. ^a). These predictions 
bracket the observed value of LKRMW and could 
even marginally explain the larger flux measured 
by IFGGT. The predicted 1^(5808) is generally 
somewhat lower than the observations. Using the 
line luminosities of Smith (1991) for WC/WO 
stars we predict W ^ 1 A; lower values are ob- 
tained using the values from the compilation of 
Schaerer &: Vacca, which accounts for differences 
between WC and WO subtypes. In view of the 
observational uncertainties (see discussion in §4) 
we conclude that the present single star models 
are able to reproduce the WR features in I Zw 18 
discussed so far with an instantaneous burst and 
a Salpeter IMF extending to Mup ^ 120-150 Mq. 

Interestingly, the predicted equivalent width 
of nebular Hell A4686 reproduces the observed 
value very well (Fig ^) . The total nebular He II 
emission in I Zw 18 is consistent with the de- 
rived number of WC stars if one assumes that 
their spectrum has the same hardness (i.e., iden- 
tical ratio of He^'''/H"'" continuum photons) as 
the well studied WO star DRl (Garnett et al. 
1991, Kingsburgh et al. 1995), and that the Ly- 
man continuum luminosity corresponds to that of 
an average WR star at low Z (Schaerer &; Vacca 
1998, their Fig. 6). Both facts support the hy- 
pothesis of Schaerer (1996) that WC/WO stars 
are responsible for nebular Hell A4686 features 
observed in extragalactic HII regions. This is in 
contrast to the result of IFGGT and LKRMW, 
who conclude, based only on the comparison of 
the observed He Il/H/3 line intensity, that the 
nebular He II emission in I Zw 18 cannot be ex- 
plained quantitatively by WC/WO stars. The 
apparent contradiction is simply explained by the 
different spatial extension of the He II and H/5 



emission seen in the HST imagery (cf. §4). The 
spectroscopic He II/H/? values in IFGGT and 
LKRMW probably overestimate the true nebu- 
lar L(He Il)/L(H/3 ) ratio because H/3 is emitted 
over a larger region than the He II. The need 
for such a geometric correction was pointed out 
by Garnett et al. (1991). An additional argu- 
ment from I Zw 18 in favor of the existence of an 
intimate link between WC/WO stars and nebu- 
lar He II emission is the spatial correlation be- 
tween the WC features and Hell A4686 found by 
LKRMW. 

5.3. WC and/or WO stars in I Zw 18 



As mentioned in §5.1 the stellar models pre- 
dict a very high C-|-0 and a low helium abun- 
dance in the WC/WO phase, which would likely 
lead to a WO instead of WC classification of 
these stars. Are WC and/or WO stars present 
in I Zw 18? A classification of the WC and WO 
subtypes observed in I Zw 18 requires more than 
the mere detection of CIV A5808 and the 4650 
bump (Smith et al. 1990; Crowther et al. 1998). 
However, as discussed by LKRMW, the observed 
FWHM of these lines is compatible with WC4 
or WC5 types (cf. Smith et al.). The width of 
CIV A5808 given by both IFGGT and LKRMW 
and the structure of the 4650 bump observed 
by the former authors are also compatible with 
those of the W03-4 star DRl (Kingsburgh et al. 
1995; Crowther et al. 1998). If the entire CIV 
A5808 emission were attributed to W03-4 stars, 
one would expect I(OV A5590)/I(CIV) « 0.15 
(Schaerer & Vacca 1998), too weak to be detected 
in the current observations. A detection of the 
broad OIV A3400, OVI A3811, or OV A5590 hues 
is required for an unambiguous detection of WO 
stars. At this time the presence of WO stars, ex- 
pected in low metallicity environments, cannot 
be excluded in I Zw 18, and the classification of 
the observed WC stars is presently not certain. 

If the fraction of evolved WR stars with (C-l-0) /He 
>1 ("WO") is low among the observed WC/WO 
stars, the present evolutionary scenario may not 
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apply and other channels leading to an earlier 
appearance of He-burning products on the stellar 
surface might be required. Two (non-exclusive) 
processes may be invoked: higher mass loss on 
the main sequence and/or the WNL phase or 
additional mixing leading to upward transport 
of processed elements. Given the uncertainties 
of the mass-loss rates of very massive stars and 
their metallicity dependence (cf. Heap et al. 1994; 
Schaerer 1998) the evolutionary scenarios of sin- 
gle stars are still uncertain. Extreme mass loss 
may also occur through Roche-lobe overflow in 
a binary system. Indeed as shown by Schaerer 
(1998) Case A binaries (Roche-lobe filled dur- 
ing MS phase) may be an important channel to 
form WR stars from high mass stars (Mjnitiai ^ 
40 Mq) in low Z environments. Rotational in- 
duced mixing can lead to an earlier entry in the 
WR phase (e.g. Meynet 1997; Meynet h Maeder 
1997). New observations of both individual mas- 
sive stars or integrated populations in low metal- 
licity environments will provide crucial informa- 
tion to guide our understanding of massive stars 
in extreme environments. 

6. Summary 

We have used an extensive collection of archival 
HST data to study the stellar content and gaseous 
distribution in I Zw 18. Continuum subtraction 
and cosmic-ray removal were found to be critical 
for the search for faint sources. Using WFPC2 
imaging and FOS spectra we were able to im- 
prove previous studies and construct a contin- 
uum free HeH A4686 map, which was used to 
search for stellar sources (WR stars) recently 
found by ground-based spectroscopy and to lo- 
cate diffuse nebular emission. The HeH A4686 
emission is associated with the NW stellar clus- 
ter, clearly displaced from the surrounding shell- 
like [O III] and Ha emission. We found evidence 
for He II sources, compatible with 5-9 WNL stars 
and/or compact nebular Hell A4686 emission, as 
well as residual diffuse emission. Only one of 
these sources is not located in NW cluster but in 



the region in between the NW and SE condensa- 
tion. 

We have done an extensive comparison be- 
tween our results and the recent ground-based 
data of IFGGT and LKRMW which allowed 
them to identify WN and WC stars in I Zw 18. 
The differences between the various data may be 
understood in terms of varying slit locations, con- 
tinuum fits and contamination by nebular lines. 
Ground-based spectra taken with better reso- 
lution and larger telescopes will allow one to 
probe the massive star content more accurately, 
to provide more stringent constraints on evolu- 
tion models, and to study the interaction be- 
tween the stars and their surrounding ISM in 
more depth. In order to avoid the difficulties in 
correlating spectra and image location, Fabry- 
Perot interferometry or other imaging spectro- 
copy would be the technique of choice. 

We have calculated evolutionary tracks for 
massive stars and synthesis models at the ap- 
propriate metallicity [Z ~ 0.02 Z©). These sin- 
gle star models predict a mass limit Mwr ~ 
9OM0 for WR stars which may become WN 
and WC/WO. For an instantaneous burst model 
with a Salpeter IMF extending up to Mup ~ 
120-150 Mq our model predictions are in rea- 
sonable agreement with the observed equivalent 
widths. The WR population in I Zw 18 has prop- 
erties consistent with those predicted by single- 
star evolution models for this metallicity. The 
WR/0 ratio forms an extension of the trend ob- 
served in WR galaxies (Maeder & Conti 1994) 
towards very low Z. Our model is also able to 
fully reproduce the observed EW of nebular Hell 
A4686 emission due to the presence of WC /WO 
stars. This quantitative agreement and the spa- 
tial correlation of nebular Hell A4686 with the 
stellar cluster and the position of WR stars fur- 
ther supports the hypothesis by Schaerer (1996) 
that WR stars are responsible for nebular He II 
emission in extra-galactic H II regions. 
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Table 1: HST archival data of I Zw 18 



Root Name 


Proposal ID 




Epoch Configuration 


Optical Element 


AA (AJ 


y2f90403t 


5434 


27 


Oct 


1994 


FOS/RD 


G190H 


722 


yii904u4t 


5434 


27 


Oct 


1994 


FOS/RD 


G190H 


722 


y2f90402t 


5434 


27 


Oct 


1994 


FOS/RD 


G570H 


2249 


u2cg0101t 


5309 


29 


Oct 


1994 


PC 


F336W 


371 


u2cg0102t 


5309 


29 


Oct 


1994 


PC 


F336W 


371 


u2cg0103t 


5309 


30 


Oct 


1994 


PC 


F336W 


371 


u2cg0301t 


5309 


30 


Oct 


1994 


PC 


"n at AWT 

F814W 


1758 


u2cg0302t 


5309 


30 


Oct 


1994 


PC 


"n a -t AWT 

F814W 


1758 


u2cg0303t 


5309 


30 


Oct 


1994 


PC 


"n at AWT 

F814W 


1758 


u2cg0401t 


5309 


31 


Oct 


1994 


PC 


F469N 


25 


u2cg0402t 


5309 


31 


Oct 


1994 


PC 


F469N 


25 


u2cgU4Uot 


5309 


31 


Oct 


1994 


PC 


r 4oyJN 


25 


u2cgU2Ult 


5309 


31 


Oct 


1994 


PC 


F555W 


1223 


u2cg0202t 


5309 


31 


Oct 


1994 


PC 


F555W 


1223 


UZCgUZUoTi 


oouy 


•^1 
ol 


Oct 


1994 


PC 


r ODD VV 


IZZO 


u2f90203t 


5434 


02 


Nov 


1994 


WF3 


F502N 


27 


u2f90204t 


5434 


02 


Nov 


1994 


WF3 


F502N 


27 


u2f90205t 


5434 


02 


Nov 


1994 


WF3 


F658N 


29 


u2f90201t 


5434 


02 


Nov 


1994 


WF3 


F702W 


1481 


u2f90202t 


5434 


02 


Nov 


1994 


WF3 


F702W 


1481 


u2f90101t 


5434 


03 


Nov 


1994 


WF3 


F702W 


1481 


u2f90102t 


5434 


03 


Nov 


1994 


WF3 


F450W 


925 


u2f90103t 


5434 


03 


Nov 


1994 


WF3 


F450W 


925 


u2f90104t 


5434 


03 


Nov 


1994 


WF3 


F555W 


1223 


u2f90105t 


5434 


03 


Nov 


1994 


WF3 


F555W 


1223 


u2cg0501t 


5309 


03 


Nov 


1994 


PC 


F656N 


22 


u2cg0502t 


5309 


03 


Nov 


1994 


PC 


F656N 


22 


u2cg0503t 


5309 


03 


Nov 


1994 


PC 


F656N 


22 


u2f90206t 


5434 


03 


Nov 


1994 


WF3 


F658N 


29 


u2f90303t 


5434 


01 


Mar 


1995 


PC 


F439W 


464 


u2f90304t 


5434 


01 


Mar 


1995 


PC 


F439W 


464 


u2f90301t 


5434 


01 


Mar 


1995 


PC 


F555W 


1223 


u2f90302t 


5434 


01 


Mar 


1995 


PC 


F555W 


1223 


u2f90305t 


5434 


01 


Mar 


1995 


PC 


F675W 


889 


u2f90306t 


5434 


01 


Mar 


1995 


PC 


F675W 


889 


y39a0305t 


6536 


08 


Jun 


1996 


FOS/RD 


G190H 


722 


y39a0306t 


6536 


08 


Jun 


1996 


FOS/RD 


G190H 


722 


y39a0303t 


6536 


08 


Jun 


1996 


FOS/RD 


G400H 


1546 


y39a0304t 


6536 


08 


Jun 


1996 


FOS/RD 


G570H 


2249 
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Table 2: Flux comparisons 



Authors^ Ha Flux [OIII] A5007 Flux 
(erg cm~^) (erg cm~^) 



This work 


3.3 


X 


10- 


-13 


2.0 


X 


10- 


-13 


DH 


4.2 


X 


10- 


-13 


2.4 


X 


10" 


-13 


DKF 


3.9 


X 


10- 


-13 


2.8 


X 


10- 


-13 


IFGGT 


1.2 


X 


10- 


-13 


8.3 


X 


10- 


-14 


PSTE 


1.8 


X 


10- 


-13 


4.6 


X 


10- 


-14 


SK'' 


1.1 


X 


nr 


-13 


7.0 


X 


10- 


-11 



°DH: Dufour & Hester (1990); DKF: Davidson et al. 
(1989); IFGGT: Izotov et al. (1997); PSTE: Pagel et al. 

(1992); SK: Skillman & Kennicutt (1993) 

'Note that only SK fluxes are reddening corrected 
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Table 3: WR stars identified in I Zw 18. 



Identification 


Position* 


Fluxes^ 


Comments 






(WNL) 




WR 


403, 465 


3.0 


in NW cluster 


WR? 


405. 464 


0.7 


in NW cluster 


WR?? 


437, 449 


0.4 


SW of the NW cluster 


WR?? 


436, 450 


0.7 


cosmic-ray contamination? 


WR?? 


435, 449 


1.2 


cosmic-ray contamination? 


WR?? 


436, 449 


2.0 


cosmic-ray contamination? 


WR 


406, 407 


1.5 


between NW and SE clusters 


1-2 WC 




LKRMW 




17±4 WNL 




IFGGT 




5±2 WC4 




IFGGT 





"x,y coordinates of each pixel in the F469N image 

''The average Une flux of 1 WNL star observed with the F469N fllter at the distance of I Zw 18 is 4.6x 10"^®erg s"^ cm"^ A"^ 



Table 4: Observed broad and nebular Hell A4686 in I Zw 18. 



Feature 






Authors'^ 




(A) 




Broad 4645 


0.029 


1.8 


LKRMW, L97 


Broad blue bump (~4686) 


0.062 




IFGGT 


Broad 5808 


0.012 - 0.031 


1.3 


LKRMW, IFGGT, L97 


Nebular Hell A4686 


0.04 


2.2 


LKRMW, IFGGT, L97 


H/3 


1.00 


56 - 127 


IFGGT, PSTE 



"LKRMW: Legrand et al. 1997, L97: Legrand (1997; private communication), IFGGT: Izotov et al. 1997, PSTE: Pagel et 
al. (1992) 
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Fig. 1.— WFPC2 V (F555W), Ha (F658N) and [OIII] A5007 (F502N) image of I Zw 18The V image 
was taken with the PC, the Ha and [OHI] images were taken with the WF3 and transformed to PC 
resolution. The rectangle delineates the area where the continuum-free helium sources were detected. 
Darker regions in the Ha image represent the strongest fluxes. [OHI] contours are from IxlO"-*^^ to 
2.6x10"^^ erg s~^ cm~^ (15 levels). Size of rectangle is 3" x 3.3" (157 pc x 172 pc). Orientation is 
given on the bottom-left corner of the V image. 
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Fig. 2. — WFPC2 V (F555W) image of I Zw 18The rectangle delineates the area in the NW region where 
the continuum-free helium sources were detected. The darkest pixels in the helium map are above the So- 
level. "WR" and "WR?" identify helium sources with fluxes equivalent to 3 WNL and 0.7 WNL stars, 
respectively. "WR??" identifies sources equivalent to 0.4, 0.7, 1.2, and 2 WNL stars. "WR" outside the 
NW region identifies a helium source equivalent to 1.5 WNL stars. 
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Fig. 3.— WFPC2 V (F555W) image of I Zw 18 
a) The white rectangles correspond to the He- 
hum sources from F469N image, b) The white 
rectangles correspond to the red stellar objects 
identified by Hunter & Thronson (1995). 



<\l' 
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Fig. 4. — Summed flux within a region of 1.5" 
along the NW region of V (F555W) and HeH 
(F469N) (pixel size is 0.046"). 
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Fig. 5. — Predicted number ratio WR 
type/(WR+0) as function of time at Z =0.0004 
for the instantaneous burst model with a Salpeter 
IMF. Sohd: total WR/(WR+0), dashed: 
WNL/(WR+0), dotted: (WC+WO)/(WR+0) 




Fig. 6. — Comparison of observed and predicted 
values WR and nebular Hell A4686 features for 
I Zw 18. Observed values (cf. Table for neb- 
ular Hell A4686 (solid triangles), broad 4645 or 
4686 emission (denoted 4650; circles), and CIV 
A5808 (squares) are shown. Differing measure- 
ments from IFGGT and LKRMW for 4650 and 
5808 are connected. Predicted values as a func- 
tion of time are given for nebular Hell A4686 
(solid line), the 4650 bump (dashed lines), and 
CIV A5808 (dotted lines) for an instantaneous 
burst model with a Salpeter IMF. The two val- 
ues for the dashed and dotted lines correspond to 
different prescriptions for the WR line luminosi- 
ties, lower panel a): equivalent widths, Upper 
panel b): line intensities relative to H/3 . 
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